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Transport and Thermodynamic Evidence for a Marginal Fermi Liquid State in ZrZn2 
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Measurements of low temperature transport and thermodynamic properties have been used to 
characterize the non- Fermi liquid state of the itinerant ferromagnet ZrZn2. We observe a T^/^ 
temperature dependence of the electrical resistivity at zero field, which becomes T^ like in an 
applied field of 9 T. In zero field we also measured the thermal conductivity, and we see a novel 
linear in T dependence of the difference between the thermal and electrical resistivities. Heat 
capacity measurements, also at zero field, reveal an upturn in the electronic contribution at low 
temperatures when the phonon term is subtracted. Taken together, we argue that these properties 
are consistent with a marginal Fermi liquid state which is predicted by a mean-field model of 
enhanced spin fluctuations on the border of ferromagnet ism in three dimensions. We compare our 
data to quantitative predictions and establish this model as a compelling theoretical framework for 
understanding ZrZn2. 

PACS numbers: 72.15.Eb, 75.40.-s, 75.50.Cc, 71.27.-l-a, 75.30.-m, 71.10.-w 



I. INTRODUCTION 



The Fermi liquid theory of the metallic state is among 
the most successful in physics, and decades of research 
have established its applicability to a wide range of sys- 
tems. The fundamental starting point of this theory is 
the existence of long-lived fermionic quasiparticles whose 
effective interactions lead only to non difFractive scatter- 
ing in the zero temperature limit and near to the Fermi 
surface. In recent years, however, there has been a prolif- 
eration of materials that display behaviour not easily un- 
derstood within the Fermi liquid picture. These so called 
'non-Fermi liquids' (NFLs) encompass a variety of sys- 
tems. Examples include the normal state of the high-Tc 
cuprates one-dimensional Tomonaga-Luttinger liq- 
uids 0, Q , paramagnetic metals near a first order phase 
transitions [1] and Kondo lattice systems @ . 

A particularly rich variety of NFL behaviour is exhib- 
ited by materials in which a magnetic phase transition is 
observed at low temperatures (see for instance reference 
0). Tuning the transition to T — i> K (a quantum crit- 
ical point) yields temperature dependences of physical 
properties that differ from the power laws predicted by 
Fermi liquid theory. Early attempts to explain this NFL 
behaviour were based on the effects of strongly enhanced 
low-frequency, long-wavelength spin fluctuations [7l-[lH 
(often called the self-consistent renormalisation or SCR 
theory). In the specific case of three-dimensional mate- 
rials near a ferromagnetic quantum critical point these 
approaches suggest characteristic transport and thermo- 
dynamic properties. For instance the electrical resistivity 
is anticipated to vary as T^^^ [13] , instead of the usual 
behaviour expected in a Fermi liquid, and the electronic 



heat capacity diverges logarithmically at low tempera- 
tures, instead of the Fermi liquid linear in temperature 
variation. 

These non-Fermi liquid temperature dependences are 
a consequence of an underlying quasiparticle scattering 
rate, r~^, which varies linearly with the excitation energy 
i? of a quasiparticle near the Fermi level. This is charac- 
teristic not of a Fermi liquid, for which varies as , 
but of a marginal Fermi liquid. The study of the marginal 
Fermi liquid (MFL) state is compelling, as it represents 
the weakest breakdown of the quasiparticle picture, and 
thus could prove to be a gateway to understanding more 
exotic departures from Fermi liquid theory. 

In this paper we examine in more detail the transport 
of heat and charge in ZrZn2 reported in a previous work 
(isj . which is a metal close to a 3D ferromagnetic criti- 
cal point. We calculate the temperature dependencies of 
transport properties using SCR theory, and show them 
to be both qualitatively and semi-quantitatively in agree- 
ment with our data. In particular we pay close attention 
to the difference between thermal and electrical resistivi- 
ties, following the analytical framework of Paglione et al. 

for clean magnetic metals. We then compare these 
results with new measurements of the electronic specific 
heat, which we argue lends further support to the validity 
of the SCR model in this material. 



II. EXPERIMENTAL DETAILS 

The intermetallic compound ZrZn2 crystallizes in the 
C15 cubic Laves structure, and is the archetypal d-band 
itinerant electron ferromagnet [T5| . Improved sample 
quality as well as reports of the possible co-existence of 
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ferromagnetism and superconductivity |16lll7| have stim- 
ulated renewed interest in this compound in recent years. 
The highest quahty samples of ZrZn2 order ferromagnet- 
ically at Tc = 28.5 K, and magnetization measurements 
reveal a small moment of 0.17 /is/Zr atom in the low- 
temperature, low-magnetic field limit which is well below 
the Curie- Weiss moment /ieS = 1.9 /is/Zr atom ^Ts']. 

The ferromagnetic transition is suppressed to zero with 
the application of only 20 kbar of hydrostatic pressure 
[13I [l8| . suggesting that ZrZn2 at ambient pressure is 
close to a ferromagnetic quantum critical point. This 
idea is supported by de Haas-van Alphen studies which 
show significantly enhanced cyclotron effective masses, as 
expected in the critical regime [Toj , and that longitudinal 
fluctuations of the local magnetisation are important [20|. 
Inelastic neutron scattering studies support the applica- 
bility of an SCR model, since the generalised magnetic 
susceptibility ixq.u) at low wavevector, q, and low fre- 
quency, Lo, has structure characteristic of overdamped 
diffusive modes with a strongly dispersive relaxation rate 
p. 

In our experiments we measured four single crystal 
samples with varying levels of impurities [H, [2^ . The 
samples were thin and platelet shaped with typical di- 
mensions of 2 mm X 0.5 mm x 0.2 mm. We label 
each crystal by its residual resistivity po and residual re- 
sistivity ratio (RRR) defined as Psook/po- These were 
measured to be po = 0.31 /zficm (RRR — 210), po = 
0.97 ^lncm (RRR = 67), po = 2.3 ^f)cm (RRR = 29) 
and po = 6.4 ^Vtcxn (RRR = 11) and were selected to 
cover a range of impurity concentrations. 

Each sample was spark cut into a convenient geometry, 
then electropolished to remove a surface layer of approx- 
imately 5 /im in depth. This step was taken to avoid sur- 
face inclusions of superconducting material, shown pre- 
viously to arise from the spark cutting procedure [17l |. 
Contacts were made in a four wire geometry using a low 
power micro-spot- welding technique [13], yielding con- 
tact resistances that were measured to be 5 mJ7 or less 
at low temperatures. 

Resistivity measurements were performed using a stan- 
dard 4-terminal low-frequency AC technique in a '^He 
cryostat and in an adiabatic demagnetisation refrigera- 
tor, with low-noise transformers to enhance the signal-to- 
noise level. Measurements were performed on the sam- 
ples with several different excitation currents between 0.1 
and 1mA and several different heating rates to ensure re- 
producibility. 

Thermal conductivity was measured with a two- 
thermometer, one-heater, steady-state technique down to 
r=0.8 K on a ^He system. The reliability of our setup 
was tested by measuring both heat and charge transport 
in silver wire, obtaining the Wiedemann- Franz law at low 
temperatures to within 2%. Measurements were checked 
using different models of Cernox thermometers and dif- 
ferent sizes of temperature gradient in order to ensure 
reliability. 

The heat capacity measurements between 0.3 and 
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FIG. 1: (color online) Main: Resistivity for the four samples 
of ZrZn2 with plotted versus T^"'^. The straight line is a fit 
to the form p(T)=po -I- AT^/^ between and 15 K. The inset 
shows the resistivity plotted against T for a wider tempera- 
ture range. 



100 K were conducted in a commercial Quantum Design 
micro calorimeter at the University of Cantabria using 
a standard relaxation technique. Only the purest sam- 
ple was measured (po = 0-31 pf2cm); a thin slab-shaped 
crystal weighing 3 mg. 



III. EXPERIMENTAL RESULTS 
A. Transport Data 

Figure [1] shows the electrical resistivity of our sam- 
ples plotted versus T^/'^ . We observe a non- Fermi liquid 
temperature dependence of the form p{T) = po + AT^^^ 
between 300 mK and approximately 15 K (or higher in 
the cleanest samples), in agreement with that observed 
in recent work [22l . |25| . No evidence of superconductiv- 
ity was observed down to 300mK in any of our crystals, 
presumably due to our chemical treatment of the surface. 

The excellent fit to a pure T^/^ power law indicates 
that the scattering of electrons by spin fluctuations is 
much larger than the scattering by phonons across a 
wide temperature range. Significant phonon scattering 
would normally be expected to produce a characteris- 
tic temperature dependence in the resistivity at low 
temperatures, crossing over to a T-linear dependence at 
higher temperatures. Neutron scattering measurements 
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FIG. 2: (color online) Electrical resistivity of ZrZn2 in a 
magnetic field. Main panel: p vs. T for the sample with 
po=0.97 ^Q.cm in ^oH = and poH = 9 T. Upper left inset: 
^.qH =9 T data, plotted vs. T^. Bottom right inset: plqH = 
T data, plotted vs. T"'^^ . The solid lines in both insets are 
fits to the data below 15 K. 



on ZrZn2 are consistent with this observation - the ob- 
served scattering intensities are up to two orders of mag- 
nitude larger than what could be expected from phonons 
alone at low temperatures (2]| . A similar conclusion was 
reached in transport studies of the cubic Laves-phase 
compounds RC02 (R=Sc,Y,Lu) by Gratz and coworkers 
[26I [27}. In these compounds, which lie near a ferromag- 
netic instability and play host to spin fluctuations, it was 
concluded that the contribution to p from phonon scat- 
tering was negligible below 30 K. 

Figure [2] shows how the T^/^ state is affected by the 
application of a magnetic field. With jii^H = 9 T, the 
resistivity recovers a characteristic Fermi liquid-like 
dependence, which may be understood to arise due to the 
suppression or gapping of the fluctuations in a magnetic 
field, as suggested by measurements of the electronic spe- 
cific heat I28| . The overall change in the resistivity with 
temperature is similarly suppressed with field, support- 
ing this interpretation. 

In can also be seen from Figure [5] that the magnetore- 
sistance of the sample is strongly temperature dependent, 
as observed previously in polycrystalline samples [l^ . At 
low temperatures, the effect is positive, as expected in an 
ordinary metal. Above T ~ 22 K however the magnetore- 
sistance becomes negative, suggesting that the decrease 
in resistivity due to the stiffening of spin fluctuations 
outweighs the increase in resistivity due to orbital niag- 
netoresistance. The recovery of power law at high 
fields helps rule out disorder scattering as the origin of 
the power law in ^lqH = 0, since this process would 
have little field dependence. 

Figure |3] shows the thermal conductivity of the three 



ZrZn2 samples. In a metallic material such as ZrZn2, the 
thermal conductivity has contributions from both elec- 
trons (kc;) and phonons (Kph)- In the ferromagnetic state 
at very low temperatures T << Tc, we might also ex- 
pect a contribution from magnons (Kmag)- A careful con- 
sideration of the relative magnitudes of each of these is 
necessary to evaluate the suitability of the MFL model 
in describing ZrZn2. 

The thermal current due to coherent, propagating, 
magnetic modes (magnons) can be difficult to distin- 
guish from other contributions to k. Such modes are 
not the overdamped, diffusive spin excitations mentioned 
previously, but coherent spin waves existing at very low 
q that may become populated in the low temperature 
regime. The contribution to thermal transport from such 
modes can be considerable, as was demonstrated in the 
antiferromagnetic insulator Nd2Cu04 [33,[3l|. However, 
in magnetic metals Kmag is limited by scattering from 
charge carriers ^32] which can greatly limit the conduc- 
tivity. We would expect Kmag to be particularly small in 
ZrZn2 due to its very small ordered moment which lim- 
its the phase space available for propagating magnons 
to g < Qsw {Qsw — — ki (X M, where k^ and k^ are 
the radii of the majority and minority spin sheets of the 
Fermi surface respectively and M is the magnetisation). 
We thus ignore the contribution from Kmag- 

Between the remaining two terms, k^i and Kph, we ar- 
gue that by far the largest contribution to thermal con- 
ductivity at low temperatures is the electronic term, as 
might be expected in metals with exceptionally low pQ. 
From Figure [31 we see that increasing the residual resis- 
tivity from 0.31 /iilcm to 6.4 fiflcm drastically reduces 
the overall thermal conductivity of the sample, and the 
characteristic peak seen in the cleanest sample at T ~ 
10 K is rapidly suppressed. Since the Kei is unlikely to 
change much between these two samples, we can infer 
that the reduction in peak height is essentially due to 
an increase in elastic scattering of electrons from impuri- 
ties. Studies of controlled levels of impurities doped into 
conventional ferromagnetic metals have shown a similar 
trend [H. 

The inset of Figure [3] shows an expanded view of the 
region near the peak observed in our cleanest sample. A 
simple linear extrapolation above and below the kink in 
the data meets at T ^ 28 K, a reasonable estimate for 
the onset temperature of the peak. This is coincident 
with the Curie temperature Tc = 28.5 K [31 , and we 
thus attribute the peak in k to the increase in electronic 
mean free path as the sample enters the ferromagnetically 
ordered state and spin fluctuation scattering is reduced. 

A more detailed treatment of the relative magnitudes 
of Kei and Kph lends quantitative support to our argu- 
ment. We can estimate Kph by comparing data from sam- 
ples with different levels of impurities as follows. We first 
define the thermal resistivity w using the Wiedemann- 
Franz law as: 
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FIG. 3: (color online) Main: Thermal conductivity of ZrZn2, 
with samples labeled by their residual resistivities. Inset: 
Zoom of the transition region in the cleanest sample, with 
po = 0.31 pfJcm. Linear fits to the data above and below Tc 
suggest a change in slope occurs at T ~28 K. 



(1) 



where the Lorentz number, Lq = l/3(7rA;B/e)^ — 2.44x 
10~^ Wri/K^. When considering the transport of both 
heat and charge by electrons, we must bear in mind that 
scattering events may affect thermal and electrical cur- 
rents in different ways. We thus define the difference 
between the thermal resistivity (w) and the electrical re- 
sistivity (p) as (5, which contains information about scat- 
tering processes: 



(2) 



The difference between the effect of scattering on the 
charge and heat current is that while charge current is 
only degraded due to the quasiparticle being deflected 
through a scattering angle 9, the heat current is also de- 
graded by any loss of energy {fuj) of the quasiparticle 
[34, ;35J. The difference between electrical and thermal 
resistivities can thus reveal information about the nature 
of inelastic scattering processes, as pointed out by Kaiser 
for the case of localized magnetic impurities in metallic 
alloys [ll] and recently used to investigate antiferromag- 
netic fluctuations in the clean magnetic metal CeRhIn5 
p^ . In the case of a ferromagnetic metal such as ZrZn2, 
such inelastic processes are attributed principally to scat- 
tering from spin fluctuations. 
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FIG. 4: The estimated phonon contribution to the total ther- 
mal conductivity using the analytical method described in the 
text. The dashed line represents a fit to Kph oc T^. The in- 
set shows the estimated phonon contribution compared to the 
total conductivity measured in the cleanest sample. 



We may now arrive at an estimate for the phonon con- 
ductivity Kpfi- At low temperatures, in keeping with 
Matthiessen's rule, we assume that the amount of im- 
purities in a given sample in the dilute limit would not 
considerably alter S. In other words, we assume that the 
elastic and inelastic scattering channels are independent. 

For two samples A and B, with differing impurity lev- 
els, we then have: 



LqT 



LoT 



Kph- 



(3) 



(4) 



We assume here that Kph will be the same in each sam- 
ple at low temperatures, since the presence of impurities 
is expected to have a minimal effect on low energy phonon 
modes [S^ . Since we measure both k and p for each sam- 
ple, at every temperature we can then solve these two 
simultaneous equations to find both S and Kph- 

The solutions of these equations for Kph are shown us- 
ing the data sets from the samples with po— 0.31 and 
6.4 pficm in Figure U) The temperature dependence of 
Kph below 25 K is reasonably close to the dependence 
typically expected from phonons scattered by electrons 
[3a |. and a fit to this form yields a coefficient Kph/T^ ~ 
0.018 mWK^'^cm^^. Using different pairs of samples in 
Equations [3] and HI yields almost identical curves for Kph, 
differing in magnitude by only 20%, which demonstrates 
the consistency of our analysis. 

The phonon contribution that we estimate using this 
method is seen to be rather small compared to the total 
conductivity of the least disorder 
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FIG. 5: The electrical and thermal resistivity of ZrZn2. The 
electrical resistivity data are shown as solid lines, while the 
thermal resistivity data are shown as data points, calculated 
using w = Lo X T / Ksi , where Ksi = /ttot - Kph as described in 
the text. 



Using our estimated form of n^h from Figure 21 we now 
subtract this from Ktot to find the electronic contribution 
to thermal transport. For the two samples with the low- 
est levels of disorder, nph is small enough so that below 
T=15 K, Ktot — Hei, meaning Kei is insensitive to the 
details of the phonon subtraction. A similar conclusion 
was drawn in thermal transport studies below 10 K in 
the very pure antiferromagnetic metal CeRhlns [l^ . 

With the above considerations we are in a position to 
study Ke/ in isolation. In Figure [5] we use k^i to calculate 
the thermal resistivity w defined in Equation [U which 
we plot alongside the electrical resistivity p. All three 
samples display the same qualitative features: w > p at 
finite temperatures, and w p in the limit T — >■ as 
expected from the Wiedemann- Franz law. 

The difference between the two resistivities is (5, de- 
fined in Equation [5] and plotted versus temperature in 
Figure [6] for the cleanest sample. The striking linear T- 
dcpendence of 5 seen in our data persists over an order 
of magnitude in temperature. In the usual metallic state 
we would expect 5{T) oc T^, since w = wq + AT"^ and 
P = Po + BT^ . The deviation from this prediction in 
ZrZn2 is clearly related to the presence of strong inelas- 
tic scattering processes operating at low temperatures. 
We rule out the possibility that 5 is influenced by the 
presence of phonon scattering at low temperatures since 
these processes lead to a very different behaviour; in the 
conventional description phonon scattering leads to a 
term in p and a term in w [H]. For T > 12 K 5{T) 
does appear to gain some upward curvature, and we take 
this temperature to be the point where phonon scatter- 
ing ceases to be negligible. We note that the same linear 
behaviour is confirmed in the sample with po=2.3 /iJlcm. 
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FIG. 6: The difference {S) of the thermal (u) = LqT / tiei) and 
electrical (p) resistivities versus temperature for the sample 
with po=0.31 /if2cm. Kej is found by subtracting nph from 
^tot , as outlined in the text. The line is a linear fit to the 
data below 12 K. 



We can compare our results with those obtained in 
a more conventional magnetic system by Paglione etal., 
who studied the antiferromagnetic metal CeRhlns with 
Tat — 3.8 K [l3l. In that study, 5 was seen to evolve 
as aT^ + bT^ below the Neel temperature and dropped 
rapidly towards zero above. The term is due to scat- 
tering from antiferromagnetic magnons. Clearly the tem- 
perature dependence of 6 at low temperatures in ZrZn2 
presents an interesting case - it is difficult to reconcile 
with either conventional electron-phonon physics or with 
the model used for CeRhlns . What we show in Section lTVl 
is that a qualitative and semi-quantitative understanding 
of 5{T) may in fact be gained by considering scattering 
from spin fluctuations on the border of ferromagnetism 
in 3D. In other words, the distinct temperature depen- 
dence of 5{T) is a characteristic of the marginal Fermi 
liquid state. 



B. Heat Capacity 

The measured specific heat capacity C(T) of ZrZn2 is 
shown in Figure [71 which are in agreement with those 
reported previously [25[ . Figure [7^ is a plot of C ver- 
sus T and Figure Wp gives C/T vs T^. The inset of 
Figure 7b gives (AC = C - Cph)/T versus T^, where 
Cph = /3T3 = {12Tr'^/5)NkB{T/eD)^ is the phonon con- 
tribution to the specific heat at low temperatures in 
which 9d is the appropriate Debye temperature. In the 
plot we have taken 9d to be 340 K, a value consistent 
with previous estimates [1^ . We show in the next section 
that the upturn of AC /T versus at low temperature is 
consistent, within experimental error, with the observed 
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FIG. 7: (a) The heat capacity, C, of ZrZn2 versus tem- 
perature, T. (b) C/T versus T^; the inset shows AC — 
{C — Cph)/T versus where Cph is a T"^ Debye-like phonon 
heat capacity whose slope is given hy do = 340 K as discussed 
in the text. 



temperature dependences of 5 and p and with the pre- 
dictions of the SCR model with parameters relevant to 
ZrZn2. 



IV. DISCUSSION 



In this section we consider the possible description of 
the transport and specific heat data in terms of the SCR 
model. We begin with an outline of the model and then 
present the predictions of the model for the model param- 
eters that we believe are relevant to ZrZn2. Finally, the 
results of the calculations are compared with our trans- 
port and specific heat data. 



A. Outline of the Self-Consistent Renormalization 
Model 

The theoretical treatment of the effects of the electron- 
electron interactions on the transport and thermody- 
namic properties of nearly or weakly ferromagnetic d- 
metals has a long history (see for instance Refs. R|, Q)- 
Here we give only an outline of the key ideas within the 
self-consistent renormalization (SCR) model mentioned 
in the introduction. 

A discussion of transport properties in this model has 
been presented most notably by Ueda and Moriya [131 ■ In 
its simplest form this model involves two isotropic bands. 
The carriers in one band, referred to as the s-band, carry 
the current while the electrons in the other band, referred 
to as c?-band, scatter the s-electrons through a local s — d 
exchange interaction. The problem can be formulated 
in terms of the scattering of s-electrons from the spon- 
taneous spin fluctuations arising in the d-electron sys- 
tem. The spin fluctuation spectrum is characterized by 
a generalized dynamical and wavevector dependent sus- 
ceptibility which emerges as a central concept in the the- 
ory. The s-electrons effectively transfer their momentum 
to the d-electrons. It is assumed implicitly that the d- 
electrons in turn ultimately transfer their momentum to 
the lattice via Umklapp processes and scattering from 
any residual impurities present. 

If we ignore residual elastic scattering due to impu- 
rities, then for T « Tp (where Tp is the Fermi tem- 
perature for the s-band), p and w in this model in the 
paramagnetic state can be written in the form 




where is the Bose function, — (e'^^/'^sT _ i 
and 

Here 77 is a measure of the s-d coupling parameter, kc 
is a wavevector cut-off, equal to the diameter of the s- 
electron Fermi surface in the usual description, and Xq,(^ 
is the generalized wavevector and frequency dependent 
magnetic susceptibility associated with the d-electrons. 
The factor (i.e. in Fp) arises from the fact that the 
current cannot be degraded by scattering of s-electrons 
through an infinitesimal angle (i.e. g 0). Importantly, 
the factor for the thermal resistivity includes a sec- 
ond, dynamical term, which allows for the degradation 
of peak current even when q is vanishingly small through 
strictly inelastic processes as discussed earlier. It is this 
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crucial term which leads to the lower temperature ex- 
ponent in w than in p. Note that if we set to 1 then 
Equation [5] is proportional to the quasiparticle scattering 
rate. 

The specific heat in this model arises primarily from 
the d-electrons, and we consider the contribution from 
dissipative modes In the paramagnetic state on the 
border of ferromagnetism at low temperature the specific 
heat may be expressed within the SCR model approxi- 
mately in the form 



c = y 

dT ^ 

q<qa 



\ nduj 



(8) 



where Si^ is the entropy of an undamped boson, i.e. an 
harmonic oscillator, which is given by 



Su = ks [(1 + nuj)\og{l + Huj) - n^logn^ 



(9) 



The parameters v and qc correspond to the degeneracy 
of spin fiuctuation modes and the wavevector cut-off re- 
spectively. In the usual description v = 3 and qc is the 
diameter of the d-electron Fermi surface. 

To evaluate the expressions for p, w and C we need 
in particular lmxq,uj, which we model using data from 
inelastic neutron scattering experiments. The cut off 
wavevectors kc and Qc can be estimated from the known 
properties of the Fermi surface, measured in quantum os- 
cillation experiments [l9l ]. The parameter rj can in prin- 
ciple be inferred from the Kadowaki- Woods relationship 
that has been shown [s^ to be consistent with the SCR 
model under conditions where this model reduces to the 
Fermi liquid limit. Since we shall be mainly interested in 
the relative magnitudes of p and w we shall for simplic- 
ity scale both quantities to a reference value of p at 15 
K (a crossover temperature above which phonon contri- 
butions to the transport properties become progressively 
less ignorable). Although we shall focus on the relative 
values of p and w, we note that by applying the above 
Kadowaki- Woods relationship we obtain in our calcula- 
tions absolute values of p and w, which are within 20% 
of the measured values for the model parameters relevant 
to ZrZn2. 

Within the SCR model Xq,i^ in the paramagnetic state 
may be expressed in the form 



X 



2 ^ 

cq ~ I — , 
79 



(10) 



where x(2^) is the static magnetic susceptibility, and c 
and 7 are constants dependent on the details of the 
d-electron band. The temperature dependence of the 
static susceptibility can in principle be determined self- 
consistently within the SCR model. Here we take x{T) 
from bulk susceptibility data and the parameters c and 
7 from inelastic neutron scattering data, which could be 
satisfactorily described by the model given by Equation 
[ini[21]. 

In the ferromagnetic state below the Curie tempera- 
ture Tc we must distinguish between the components of 



the spin fluctuations parallel (||) to the ordered moment 
and perpendicular (J-) to the ordered moment. The prin- 
cipal changes are as follows. Firstly, the static suscepti- 
bility x{T) becomes strongly dependent on the compo- 
nent (there are two transverse (_L) components and one 
longitudinal (||) component). Secondly, the component 
of the generalized susceptibility has a spin- wave form be- 
low the cut-off qsw defined earlier. In this regime we have 



where 



xll^ = cq\i-^^/nl) 



Mlq = Dq^ = 2^iBMcq^ 



(11) 



(12) 



is the spin wave spectrum at small q. 

From inelastic scattering data for ZrZn2 we have ap- 
proximately 7=1.8 peVk and c= 3x10^ |21|]. From 
the measured Fermi surface of ZrZn2 we estimate qsw — 
0.07 [l^. From bulk measurements of the mag- 

netic equation of state we arrive at models for x±{T) 
and X\\iT) as shown in Figure H] and M{T) = A/o(l - 
rV^c)^^^ where Mq =31 G. We note that the resuffs 
of the calculations are not sensitive to precise forms of x 
and M or even to the precise values of the other param- 
eters. For non-Fermi liquid behaviours to arise in this 
model we require principally that x~^ is small compared 
with cg^ for the characteristic wavevectors for thermally 
excited spin fiuctuations. In other words, we require that 
the characteristic thermal wavevector, qT, is large com- 
pared with the magnetic correlation wavevector (i.e., the 
inverse of the magnetic correlation length ^ = \/xfc) [Ij . 

In practice, d-metals such as ZrZn2 cannot be repre- 
sented simply in terms of two distinct species of electron 
as the above model suggests and the terms "s-electron" 
and "d-electron" may not have precise meanings. It is of 
interest, however, to see if the model can capture some of 
the behaviours of p, w and C, which we have measured 
in ZrZn2. 

It has been emphasized in more recent theoretical 
works that there should be non-anal ytic corrections to 
certain aspects of the SCR model [39l - l4ll |. Although rel- 
atively weak in 3D, they can nevertheless strongly modify 
the magnetic equation of state in particular and give rise 
to a first-order ferromagnetic transition near the quan- 
tum critical point and potentially to a magnetically in- 
homogeneous state. We shall therefore not attempt to 
calculate the magnetic equation of state self-consistently 
in this paper and take the temperature dependence of as- 
pects of the magnetic equation of state from experiment. 



B. The Marginal Fermi Liquid Limit of the SCR 
Model 

In the limit qr » it may easily be shown from 
Equations [5] -[TOl that in 3D 



p oc 



rp5/3 



(13) 
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FIG. 8: (color online) Plots showing the form of the di- 
mensionless inverse volume susceptibility (in c.g.s units of 
Oe/ (emu cm"'')) used in the resistivity calculations for ZrZn2. 



The expressions of these are as follows: X|[^ = 2(36rc(l - 
(T/Tc)^/^)) and xl' = for T < Tc and x^' = Xl' = 
70(r - Tc) for T > Tc- Here we take Tc =28.5 K. 

S = (w ~ p) <xT (14) 

oc T\og{T*/T) (15) 

and 

C cx T\og{T*/T) (16) 



where T* is an appropriate temperature scale in each case 
separately. Note that the temperature dependence of the 
scattering rate is different from that of 5 due to the 
{lo/TY factor in Equation [T] Thus, in this limit, the 
SCR model reduces to the MFL model discussed in the 
introduction, while in the opposite limit, qT << the 
SCR model reduces to the FL model. For the parameters 
relevant to ZrZn2 we find that the crossover to the FL 
model occurs only at temperatures below 1 K at ambient 
pressure [i^l. 

The characteristic marginal Fermi liquid exponents for 
p and S can be understood by analogy to the correspond- 
ing electron-phonon scattering problem. In that problem 
p oc q!^~^^ and 6 (x q!^ , where d is the dimension of space. 
The extra factor of q^ in p comes from the q^ term in Fp 
(Equation |6]) which is absent in the difference between w 
and d. The temperature dependence of qx is given by the 
dynamical exponent z via T on q^ so that one expects 

ptxT— , (17) 
S(xTi. (18) 



For the electron-phonon scattering problem z = 1 and 
in 3D we recover the usual results, p (x and 6 oc T^. 
For electron paramagnon scattering, however, z = 3 (see 
EquationlTUlin the limit 0) and so in 3D we expect 

p oc and (5 oc T as in Equations [13] and [H We 

caution that for dissipative modes such as paramagnons 
this elementary treatment fails if z is unity or lower. In 
this case, p and S are determined by spin-fluctuations 
over all of g-space rather than just at low q and to obtain 
the FL limit of the SCR model one must return to the 
full expressions for p and w in Equation ([S]). 



C. Calculations of the Electrical and Thermal 
Resistivities 

Using the model and parameters given above we can 
calculate the electrical and thermal resistivities and the 
difference 5 = w — p. The results of these calculations 
are shown in Figure [3] for several values of fee (which in 
the model represents the typical diameter of the Fermi 
surface of the light electrons that principally carry the 
current). 

The value of fee that gives the correct magnitude of 
the scaled quantity 6{T)/ p{lb K) turns out to be around 
fcc=0.4 A~^. It is interesting to note that this is close to 
the observed diameter of the sheet of the Fermi surface 
of ZrZn2 with the highest characteristic Fermi velocity 
observed in this system. However, this fitted value of fee 
should not be taken too seriously, given the complexity 
of the Fermi surface of ZrZn2 and the relative simplic- 
ity of the present model. What seems to be significant, 
however, is that the calculated temperature dependence 
of p and 5 are qualitatively very similar to the observed 
values given in Section IIII Al irrespective of the precise 
value chosen for fcc, or indeed of the precise values of the 
other model parameters. 

Within the experimental temperature range of inter- 
est, we find that the SCR model reduces essentially to the 
MFL model as discussed in the introduction. In partic- 
ular, p is predicted to vary approximately as T^/'^ and 5 
approximately as T in the temperature range below 15 K 
(above which phonon contributions become progressively 
more important in the measurements) and above about 
1 K (well below which the SCR model predicts a Fermi 
liquid form for p and w). 



D. Calculations of the Specific Heat 

Using the model and parameters in Section IIV Al we 
can also calculate the specific heat contribution of the 
spin-fluctuations. Below Tc we consider the contribution 
from both dissipative modes and from spin waves with 
q < <lsw The results of the calculations are shown in 
Figure [TOl in the form of C/T vs for comparisons with 
the inset of Figure[7lD of the measured specific heat minus 
the estimated phonon contribution. The calculated heat 
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FIG. 9: (color online) Plot showing (a) the calculated electri- 
cal resistivity plotted against t'^/^ and (b) the calculated dif- 
ference between the thermal and electrical resistivities plotted 
against T. The results are normalised relative to the calcu- 
lated electrical resistivity at 16 K. The results are shown for 
three values of the wavevector cut-off parameter (0.4 
(bottom), 0.6 A"^ (middle) and 0.8 A"^ (top)). The data for 
kc = 0.6 A~^ and fee = 0.8 A~^ are shifted upwards by 0.1 
and 0.2 respectively, for clarity. 

capacity depends approximately logarithmically on the 
upper cut-off qc and is also somewhat sensitive to the 
spin- wave cut-ofF qsw , although the relative contribution 
from the dissipative modes with q > qsw dominates at all 
temperatures. The results in Figure [10] correspond to the 
choice qc— 1 A~^, which is approximately equal to the 
observed diameter of the large sheets of the Fermi surface 
with high effective mass (and so low Fermi velocity). 

To check on the sensitivity of the result to the lower 
spin wave cut-ofF several values of qsw around the value 
of 0.07 A-i given in Section llV Al have been used. 

Comparing the inset of Figure. Wp with Figure [TOl we 
see that the measurements are qualitatively consistent 
with the predictions of the SCR model, especially for a 
reasonable value for qsw In particular, the upturn in 
C /T with decreasing temperature in the inset of Figure 
7b is consistent with the predictions of the SCR model. 
In the limit qT » (see Section liVB[) this upturn is 
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FIG. 10: Plot showing the calculated spin fluctuation heat 
capacity for several values of qaw, namely 0.05 A~^, 0.08 A~^ 
and 0.1 A"^ The effect of increasing qsw is to cut off the low 
temperature divergence of C/T. 

reminiscent of the logarithmic term in C expected for a 
marginal Fermi liquid. The magnitude of the calculated 
C /T is somewhat lower than the measured value. This is 
reasonable since the SCR model does not include the ef- 
fects of the electron-phonon interaction, nor, in the form 
defined in Section flV Al the effects of any antifcrromag- 
netic spin fluctuations that may be present. 



V. CONCLUSIONS 

In summary, we find that the temperature dependences 
of the electrical and thermal resistivities and the elec- 
tronic specific heat in ZrZn2 can be understood qualita- 
tively and even semi-quantitatively in terms of the self- 
consistent renormalization model with model parameters 
inferred from independent measurements of the magnetic 
equation of state, the inelastic neutron scattering cross- 
section and the Fermi surface. Over a wide range in tem- 
perature, the SCR model reduces to the marginal Fermi 
liquid model characterized by an electrical resistivity p 
varying as T^/'^, a thermal resistivity w - p varying as 
T and a logarithmic divergence in the low temperature 
specific heat. We believe that taken together, these find- 
ings provide the most compelling evidence thus far for 
the existence of the marginal Fermi liquid state in the 
field of itinerant-electron ferromagnetism. 
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